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A SIMPLE ELEMENT FOR MULTILAYER BEAMS 
IN NASTRAN THERMAL STRESS ANALYSIS 

W.T. Chen and S.K. Wadhwa 
International Business Machines Corporation 

I. INTRODUCTION 


In the study of structural integrity for electronic packages under thermal 
cycles, an important consideration is the distortion due to the differential 
thermal expansions of multilayer dissimilar materials. A simple antecedent for 
such cases Is the classical bimetallic beam analysis by S, Timoshenko. In the 
application of NASTRAN, such structural members are usually represented by 
bar elements with multi-point constraint cards to enforce the interface conditions. 
While this is a very powerful method in principle, one finds that in practice 
the process for specification of constraints becomes tedious and error prone, 
unless the geometry is simple and the number of grid points low. An alterna- 
tive approach has been found within the framework of the NASTRAN program. 
This approach makes use of the idea that a thermal distortion in a multilayer 
beam may be similar to a homogenous beam with a thermal gradient across the: 
cross section. This paper contains the exact mathematical derivation for the 
equivalent beam, and all the necessary formulae for the equivalent parameters 
in NASTRAN analysis. Some numerical examples illustrate the simplicity and 
ease of this approach for finite element analysis such as NASTRAN. 

II . ANALYSIS 

Consider an n-layer composite beam of dissimilar materials at constant 
temperature T, having an external loading N as the axial force, and M as 
the bending moment. The cross section is s?iown in Figure 1. ^ 



tj=: Thickness of layer 
Zj= Coordinate of interface 

Fig. 1 - Cross Section of n layer Composite Beam 
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ti, t 2 ..t are the thicknesses of various layers and are the coordinates 

at the interface. With the coordinates as specified, ' 



Strain at any point can be written as 


z-d 

e = e + — ~ — “ hT 

XX t" 


where e = elastic strain 

XX 

e = average axial strain of composite 
z = coordinate of the point 

d = composites neutral axis 

r = radius of curvature 

aT = thermal strain 

a = thermal coefficient of expansion of the layer 



Let Ej, a., t. be the Young's modulus, thermal coefficient of expansion 
and thickness of ilh layer with uniform width b. 

Define z =0 and/ 
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then 
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be shown as 
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From Equation (3) , since N should be independent of r, 
therefore 

pd - y = o 

d = — = distance of neutral axis. 


substituting this in Equations (8) and (9) we get 



Thus the external force and moment have been derived for n layer structure. 
Expression 6-y* is simply area stiffness 0n which E is included^ moment of 

J 

inertia about centroid, and can be proved as follows. 

6 = Area Stiffness M.i. about coordinaie axIsQas defined by Equation 5^ . 
Ml about centroid fusing parallel axis theorem^ = 


,6 

6 


= 6 


- Area Stiffness 

- P . y2 
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X d^ 
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Consider now an equivalent homogenous beam of area A, centroid moment of 
inertia 1, Young's modulus E, thermal coefficient of expansion having thermal 
gradient dT . 
dz 

Then 


G 


XX 


e + 



aT 


where the coordinate axes are the neutral axes of the beam. By using the 
previous formulation it can be derived that 

= AEe - AEaT 

(12) 

M = J_ - la dT 

y r — 

dz 

(13) 

Comparing Equation 10 to 12 and 11 to 
ships can be written for equivalency. 

13, the following parametric relation 

AE = [1 

(14) 

AEaT = 

(15) 

1=6 ^ 

P 

la N 

dz T Q T 

(16) 

(17) 


Thus assuming E for the equivalent beam. Area A, a , I, dT can be calculated 
from Equations 14 through 17. dz 


Thus a nonhomogenous beam [[specified as a finite element in NASTRAl^ 
at constant temperature can be solved by a mathematical equivalent homogenous 
beam with temperature gradient at the cross section. 

III. TIMOSHENKO'S BIMETALLIC BEAM 

It can further be proved that for two layer structure, in the absence of 
external forces, Equation 9 degenerates to Timoshenko bimetallic beam, 

Consider now a bimetallic beam of thickness ti= tz= t, and therefore 
2l = t, Z2 = 2t, 
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From Equations 3, 4, 5, 6 and 7 


|i ;= (El 4- Ea) 

2 

Y = — [Ei + 3Ez] 

I 

8 = 4- 7eD 

N-r = bT y- (Eitti + EzCk) 

(Eitti + 3E2ttz) 


M 


T 



T 



T 






ExEa 

Ei+Ez 



[E. + 7E.]- ^ 


(Ei + 3E2)^ 
Ei+ E 2 



(E 1 +E 2 ) + 

O 



Substituting in Equation 11, we get 




(Ei+Ei) + 



' El Eg 

.E 1 +E 2 



{03 - ai ) 


El Eg 
Ei+ Eg 


bt^ ETEg 
~ E 1 +E 2 



(02 - Oi) T 


In the absence of external force, the expression in parentheses Is zero, which 
is identical to Timoshenko's bimetallic beam. 
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IV. EXAMPLE 


For the purpose of illustration, the example in Strength of Materials 
Part I by TimoshenkojjJ will be discussed. To model a bimetallic beam 
section, one defines four grid points and two bar elements. For each pair of 
grid points, continuity of the interface displacements must be specified. 
Multipoint constraints or MPC cards must be used to relate the displacements 
and rotation of one bar element to the other of the following form 

Ui + hi 01 = U2 + hj 02 

It should be noted that this could become quite complicated when many elements 
are used to represent a curved beam '.n space. 

An equivalent beam with the same mechanical response will now be 
constructed. The bending curvature of the bimetallic beam will now be induced 
by the TEMPRB Bulk Data card. The bending stress in the original bimetallic 
beam due to differential expansion will now be replaced by the thermal moment 
thermal gradient defined by the TEMPRB card for each bar element. In the 
example in reference , page 219-, we shall use 

b = h = 1, El = 1 X 10\ C 2 = 1.15 X 1Q\ Oi = 1 x lO"’ , = 2 x 10"‘* 

Then it turns out that the equivalent beam should be 

A = 1.8695, ! = .62091, E = 1.15 x 107, a = 1.5348 x 10"\- 
with 100 degree rise in temperature, the gradient should be 48.8. 


V. DISCUSSION 

The method proposed in this paper has been found to be very useful in 
analysis dealing with deformation associated with multi-layered curved beam 
structures undergoing thermal loads. The main advantages are (!) elimination 
of the time-consuming task in specifying multipoint constraints, and (2) 
reduction in number of grid points. 
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